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Abstract 

Recent development of trilayer graphene nanoribbon Schottky-barrier field-effect transistors (FETs) will be governed 
by transistor electrostatics and quantum effects that impose scaling limits like those of Si metal-oxide-semiconductor 
field-effect transistors. The current-voltage characteristic of a Schottky-barrier FET has been studied as a function of 
physical parameters such as effective mass, graphene nanoribbon length, gate insulator thickness, and electrical 
parameters such as Schottky barrier height and applied bias voltage. In this paper, the scaling behaviors of a 
Schottky-barrier FET using trilayer graphene nanoribbon are studied and analytically modeled. A novel analytical 
method is also presented for describing a switch in a Schottky-contact double-gate trilayer graphene nanoribbon 
FET. In the proposed model, different stacking arrangements of trilayer graphene nanoribbon are assumed as metal 
and semiconductor contacts to form a Schottky transistor. Based on this assumption, an analytical model and 
numerical solution of the junction current-voltage are presented in which the applied bias voltage and channel 
length dependence characteristics are highlighted. The model is then compared with other types of transistors. The 
developed model can assist in comprehending experiments involving graphene nanoribbon Schottky-barrier FETs. 
It is demonstrated that the proposed structure exhibits negligible short-channel effects, an improved on-current, 
realistic threshold voltage, and opposite subthreshold slope and meets the International Technology Roadmap for 
Semiconductors near-term guidelines. Finally, the results showed that there is a fast transient between on-off states. 
In other words, the suggested model can be used as a high-speed switch where the value of subthreshold slope is 
small and thus leads to less power consumption. 

Keywords: Trilayer graphene nanoribbon (TGN), ABA and ABC stacking, TGN Schottky-barrier FET, High-speed 
switch 



Background 

Graphene, as a single layer of carbon atoms with hex- 
agonal symmetry and different types such as monolayer, 
bilayer, trilayer, and multilayers, has attracted new 
research attention. Very high carrier mobility can be 
achieved from graphene-based materials which makes 
them a promising candidate for nanoelectronic devices 
[1,2]. Recently, electron and hole mobilities of a sus- 
pended graphene have reached as high as 2 x 10 5 cm 2 / 
V-s [3]. Also, ballistic transport has been observed at 
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room temperature in these materials [3]. Layers of gra- 
phene can be stacked differently depending on the hori- 
zontal shift of graphene planes [4,5]. Every individual 
multilayer graphene sequence behaves like a new mater- 
ial, and different stacking of graphene sheet lead to dif- 
ferent electronic properties [3,6,7]. In addition, the 
configuration of graphene layers plays a significant role 
to realize either metallic or semiconducting electronic 
behavior [4,8,9]. 

Trilayer graphene nanoribbon (TGN), as a one- 
dimensional (ID) material, is the focus of this study. The 
quantum confinement effect will be assumed in two 
directions. In other words, only one Cartesian direction 
is greater than the de Broglie wavelength (10 nm). As 
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shown in Figure la, because of the quantum confine- 
ment effect, a digital energy is taken in the y and z direc- 
tions, while an analog type in the x direction. It is also 
remarkable that the electrical property of TGN is a 
strong function of interlayer stacking sequences [10]. 
Two well-known forms of TGN with different stacking 
manners are understood as ABA (Bernal) and ABC 
(rhombohedral) [11]. The simplest crystallographic 
structure is hexagonal or AA stacking, where each layer 
is placed directly on top of another; however, it is un- 
stable. AB (Bernal) stacking is the distinct stacking struc- 
ture for bilayers. For trilayers, it can be formed as either 
ABA, as shown in Figure 1, or ABC (rhombohedral) 
stacking [1,12]. Bernal stacking (ABA) is a common hex- 
agonal structure which has been found in graphite. 
However, some parts of graphite can also have a 
rhombohedral structure (the ABC stacking) [6,13]. The 
band structure of ABA-stacked TGNs can be assumed 
as a hybrid of monolayer and bilayer graphene band 
structures. The perpendicular external applied electric 
or magnetic fields are expected to induce band crossing 
variation in Bernal-stacked TGNs [14-16]. Figure 1 indi- 
cates that the graphene plane being a two-dimensional 
(2D) honeycomb lattice is the origin of the stacking 
order in multilayer graphene with A and B and two non- 
equivalent sublattices. 

As shown in Figure 1, a TGN with ABA stacking has 
been modeled in the form of three honeycomb lattices 
with pairs of equivalent sites as {A^BJ, {A 2 ,B 2 }, and {A 3 , 
B 3 } which are located in the top, center, and bottom 
layers, respectively [11]. An effective-mass model utiliz- 
ing the Slonczewski- Weiss -McClure parameterization 
[17] has been adopted, where every parameter can be 
compared with a relevant parameter in the tight-binding 
model. The stacking order is related to the electronic 
low-energy structure of 3D graphite-based materials 
[18,19]. Interlayer coupling has been found to also affect 



the device performance, which can be decreased as a re- 
sult of mismatching the A-B stacking of the graphene 
layers or rising the interlayer distance. A weaker inter- 
layer coupling may lead to reduced energy spacing be- 
tween the subbands and increased availability of more 
subbands for transfer in the low-energy array. Graphene 
nanoribbon (GNR) has been incorporated in different 
nanoscale devices such as interconnects, electromechan- 
ical switches, Schottky diodes, tunnel transistors, and 
field-effect transistors (FETs) [20-24]. The characteristics 
of the electron and hole energy spectra in graphene cre- 
ate unique features of graphene-based Schottky transis- 
tors. Recently, the fabrication and experimental studies 
as well as a hypothetical model of G -Schottky transistors 
have been presented [25]. The studies have focused to- 
wards the properties of TGN, and a tunable three-layer 
graphene single-electron transistor was experimentally 
realized [6,26]. 

In this paper, a model for TGN Schottky-barrier (SB) 
FET is analyzed which can be assumed as a ID device 
with width and thickness less than the de Broglie wave- 
length. The presented analytical model involves a range 
of nanoribbons placed between a highly conducting sub- 
strate with the back gate and the top gate controlling 
the source-drain current. The Schottky barrier is defined 
as an electron or hole barrier which is caused by an elec- 
tric dipole charge distribution related to the contact and 
difference created between a metal and semiconductor 
under an equilibrium condition. The barrier is found to 
be very abrupt at the top of the metal due to the charge 
being mostly on the surface [27-31]. TGN with different 
stacking sequences (ABA and ABC) indicates different 
electrical properties, which can be used in the SB struc- 
ture. This means that by engineering the stack of TGN, 
Schottky contacts can be designed, as shown in Figure 2. 
Between two different arrangements of TGN, the semi- 
conducting behavior of the ABA stacking structure has 
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Figure 1 TGN. (a) As a one-dimensional material with quantum confinement effect on two Cartesian directions, (b) ABA-stacked [17]. 
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turned it into a useful and competent channel material 
to be used in Schottky transistors [32]. 

In fact, the TGN with ABC stacking shows a semi- 
metallic behavior, while the ABA-stacked TGN shows a 
semiconducting property [32]. A schematic view of TGN 
SB FET is illustrated in Figure 3, in which ABA-stacked 
TGN forms the channel between the source and drain 
contacts. The contact size has a smaller effect on the 
double-gate (DG) GNR FET compared to the single-gate 
(SG) FET. 

Due to the fact that the GNR channel is sandwiched 
or wrapped through by the gate, the field lines from the 
source and drain contacts were seen to be properly 
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Figure 3 Schematic representation of TGN SB FET. 



screened by the gate electrodes, and therefore, the 
source and drain contact geometry has a lower impact. 
The operation of TGN SB FET is followed by the cre- 
ation of the lateral semimetal-semiconductor-semimetal 
junction under the controlling top gate and relevant en- 
ergy barrier. 

Methods 

TGN SB FET model 

The scaling behaviors of TGN SB FET are studied by 
self-consistently solving the energy band structure equa- 
tion in an atomistic basis set. In order to calculate the 
energy band structure of ABA-stacked TGN, the 
spectrum of full tight-binding Hamiltonian technique 
has been adopted [33-37]. The presence of electrostatic 
fields breaks the symmetry between the three layers. 
Using perturbation theory [38] in the limit of v F \k\ « V« 
t ± gives the electronic band structure of TGN as [35,39] 



E k = ±(ak-/3k 3 ), 



(1) 



where k is the wave vector in the x direction, a - 



Vf.V 



iW2'P = iWw' tjL is the ho PP in g ener gy> v f is the Fermi 
velocity, and V is the applied voltage. The response of 
ABA-stacked TGN to an external electric field is differ- 
ent from that of mono- or bilayer graphene. Rather than 
opening a gap in bilayer graphene, this tuned the magni- 
tude of overlap in TGN. Based on the energy dispersion 
of biased TGN, wave vector relation with the energy (£- 
k relation) shows overlap between the conduction and 
valence band structures, which can be controlled by a 
perpendicular external electric field [6,39]. The band 
overlap increases with increasing external electric field 
which is independent of the electric field polarity. 
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Figure 4 The DOS of the TGN with ABA stacking. 



Moreover, it is shown that the effective mass remains 
constant when the external electric field is increased 
[3,33]. As an essential parameter of TGNs, density of 
states (DOS) reveals the availability of energy states, 
which is defined as in [40,41], To obtain this amount, 
derivation of energy over the wave vector is required. 
Since DOS shows the number of available states at each 
energy level which can be occupied, therefore, DOS, as a 
function of wave vector, can be modeled as [39] 



DOS(£) 



c(de + Vf + e 2 ) 



(2) 



where E is the energy band structure and A, B, C, D, 
and F are defined as A = -6.2832a, B = 14.3849a 2 ^, 
C = ™, D = -91?, and F = =Mj9Qc?. As shown in 
Figure 4, the DOS for ABA-stacked TGN at room 
temperature is plotted. As illustrated, the low-DOS 
spectrum exposes two prominent peaks around the 
Fermi energy [39]. 

The electron concentration is calculated by integrating 
the Fermi probability distribution function over the 
energy as in [42]. Biased ABA-stacked TGN carrier 
concentration is modified as [43] 



where x = the normalized Fermi energy is rj = E ^ + j f 

F 



k B T 

and M and N are defined as M = and N 

Kb J- \KbT) 

Based on this model, ABA-stacked TGN carrier con- 
centration is a function of normalized Fermi energy 
{rj). The conductance of graphene at the Dirac point 
indicates minimum conductance at a charge neutrality 
point which depends on temperature. For a ID TGN FET, 
the GNR channel is assumed to be ballistic. The current 
from source to drain can be given by the Boltzmann trans- 
port equation in which the Landauer formula has been 
adopted [44,45]. The number of modes in corporation with 
the Landauer formula indicates conductance of TGN that 
can be written as [32] 




(4) 



where the momentum (k) can be derived by using Carda- 
no's solution for cubic equations [46]. Equation 4 can be 
assumed in the form G = NiGi + N 2 G 2 , where Ni = 2aq 2 / 
Ih and A/2 = -6^/lh. Since Gi is an odd function, its value 
is equivalent to zero. Therefore, G = N 2 G 2 [32], where 




+ 




aiL \\ + e k B T j 



\dE. 



(5) 



This equation can be numerically solved by employing the 
partial integration method and using the simplification 
form, where x = (E - A)//c B rand rj = (£ F - A)/k B T. Thus, 
the general conductance model of TGN will be obtained 
[32] as 



(ksT)dx 



(k B T)i \p{x+M) + ^N+{x+M) 2 } 3 



C(k B Tf { D(x +M) + y N + (x + M) 




(3) 



exp(# — rj)) 
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It can be seen that the conductivity of TGN increases 
by raising the magnitude of gate voltage. In the Schottky 
contact, electrons can be injected directly from the metal 
into the empty space in the semiconductor. When elec- 
trons flow from the valence band of the semiconductor 
into the metal, there would be a result similar to that for 
holes injected into the semiconductor. So, the establish- 
ment of an excess minority carrier hole in the vicinity is 
observed [28]. The current moves mainly from the drain 
to the source which consists of both drift and diffusion 
currents. The created 2D anticipated framework is 
expected to cause an explicit analytical current equation 
in the subthreshold system. Considering the weak inver- 
sion region, the diffusion current is mainly dominated 
and relative to the electron absorption at the virtual 
cathode [47]. A GNR FET is a voltage-controlled tunnel 
barrier device for both the Schottky and doped contacts. 

The drain current through the barrier consists of ther- 
mal and tunneling components [48]. The effect of 
quantum tunneling and electrostatic short channel is 
not treated, which makes it difficult to study scaling 
behaviors and ultimate scaling limits of GNR SB FET 
where the tunneling effect cannot be ignored [20]. The 
tunneling current is the main component of the whole 
current which requires the use of the quantum trans- 
port. Close to the source within the band gap, carriers 
are injected into the channel from the source [49]. In 
fact, the tunneling current plays a very important role in 
a Schottky contact device. 

The proposed model includes tunneling current 
through the SB at the contact interfaces, appropriately 
capturing the impact of arbitrary electrical and physical 
factors. The behavior of the proposed transistor over the 
threshold region is obtained by modulating the tunnel- 
ing current through the SBs at the two ends of the chan- 
nel [20]. The effect of charges close to the source for a 
SB FET is more severe because they have a significant 
effect on the SB and the tunneling possibility. When the 
charge impurity is situated at the center of the channel 
of a SB FET, the electrons are trapped by the positive 
charge and the source-drain current is decreased. If the 



charges are situated close to the drain, the electrons will 
collect near the drain. In this situation, low charge dens- 
ity near the source decreases the potential barrier at the 
beginning of the channel, which opens up the energy 
gap more for the flow of electrons from the source to 
the channel [50]. 

Electrons moving from the metal into the semicon- 
ductor can be defined by the electron current density 
/ m ^ s , whereas the electron current density / s ^ m refers to 
the movement of electrons from the semiconductor into 
the metal. What determines the direction of electron 
flow depends on the subscripts of the current. In other 
words, the conventional current direction is opposite to 
the electron flow. / s ^ m is related to the concentration of 
electrons with velocity in the x direction to subdue the 
barrier [28]: 

Js^m = e j v x dn, (7) 

— oo 

where e is the magnitude of the electronic charge and v x 
is the carrier velocity in the direction of transport: 

/total — Jm^s ~, • (8) 



High carrier mobility reported from experiments on 
graphene leads to assume a complete ballistic carrier 
transport in the TGN, which means that the average 
probability of injected electron at one end that will 
transmit to the other end is approximately equal to 1: 

/total = Jm^s = Js^m • (9) 

Kinetic energy, as a main parameter, is considered over 
the Fermi level, and the current density-voltage response 
of the TGN SB FET device is determined with respect to 
the carrier density and its kinetic energy as 
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(10) 



where rj~ ^ r ^ r (Va is the applied bias voltage and V T is 
the thermal voltage) [51]. The dependence of the drain 
current on the drain-source voltage is associated with 
the dependence of rj on this voltage given by 



(^gt - V{y))e 
KrT 



(ii) 



where V GT = V GS - V T and V(y) is the voltage of chan- 
nel in the y direction. By solving Equation 11, the nor- 
malized Fermi energy can be defined as 



K B T 



Vgt Vbs 



V 2 ' 



(12) 



In order to obtain an analytical relation for the contact 
current, an explicit analytical equation for the electric po- 
tential distribution along the TGN is presented. The chan- 
nel current is analytically derived as a function of various 
physical and electrical parameters of the device including 
effective mass, length, temperature, and applied bias volt- 
age. According to the relationship between a current and 
its density, the current-voltage response of a TGN SB 
FET, as a main characteristic, is modeled as 



Vlel f 
\[m* J 



(k B Tfxidx 



(k B T)l { D{x+M) + ^N+{x+M) 2 } 3 



C{k B Tf^D{x + M) + yJN + 0 + M) 2 
l + exp(*-^[y GT y DS -Y])) 



(13) 



where / is the length of the channel. 



Results and discussion 

In this section, the performance of the Schottky-contact 
double-gate TGN FET is studied. A novel analytical 
method is introduced to achieve a better understanding of 
the TGN SB switch devices. The results will be applied to 
identify how various device geometries provide different 
degrees of controlling transient between on-off states. The 



numerical solution of the presented analytical model in 
the preceding section was employed, and rectification 
current-voltage characteristic of TGN SB FET is plotted 
as shown in Figure 5. 

It further revealed that the engineering of SB height does 
not alter the qualitative ambipolar feature of the current- 
voltage characteristic whenever the gate oxide is thin. The 
reason is that the gate electrode could perfectly screen the 
field from the drain and source for a thin gate oxide (less 
than 10 nm). The SB whose thickness is almost the same 
as the gate insulator diameter is almost transparent. How- 
ever, the ambipolar current-voltage (/-V) characteristic 
cannot be concealed by engineering the SB height when 
the gate insulator is thin. Lowering the gate insulator thick- 
ness and the contact size leads to thinner SBs and also 
greater on-current. Since the SB height is half of the band 
gap, the minimum currents exist at the gate voltage of 
^G,min = 1/2 Vd, at which the conduction band that bends 
at the source extreme of the channel is symmetric to the 
valence band and also bends at the drain end of the chan- 
nel, while the electron current is the same as the hole 
current. The consequence of attaining the least leakage 
current is the same as TGN SB FET with middle-gap SBs 
[23]. Raising the drain voltage leads to an exponential in- 
crease of the minimal leakage current which shows the im- 
portance of proper designing of the power supply voltage 
to ensure small leakage current. As depicted in Figure 6, 




0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
V DS (Voltage) 

Figure 5 Simulated / D (uA) versus V DS (V) plots of TGN 
Schottky-barrier FET (L = 25 nm, V GS = 0.5 V). 
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Figure 6 / D (uA)-l/ DS (V) characteristic of TGN SB FET at different values of V GS for L = 



100 nm. 



the proposed model points out strong gate-source voltage 
dependence of the current-voltage characteristic showing 
that the Vgs increment effect will influence the drain 
current In other words, as V GS increases, a greater value 
of 7 D results. As the drain voltage rises, the voltage drop 
through the oxide close to the drain terminal reduces, and 
this shows that the induced inversion charge density close 
to the drain also decreases [28] . The slope of the 7 D versus 
Vds curve will reduce as a result of the decrease in the in- 
cremental conductance of the channel at the drain. This 
impact is indicated in the I^-Vds curve in Figure 6. If V DS 
increases to the point that the potential drop across the 
oxide at the drain terminal is equal to V T , the induced 



inversion charge density is zero at the drain terminal. At 
that point, the incremental conductance at the drain is nil, 
meaning that the slope of the IwVds curve is zero. We 
can write 



F DS (sat) = V GS ~ V T , 



(14) 



where Vd S (sat) is the drain-to-source voltage which is cre- 
ating zero inversion charge density at the drain terminal. 
When V DS is more than the V DS (sat) value, the point in 
the channel where the inversion charge is zero moves 
closer to the source terminal [28]. In this case, electrons 
move into the channel at the source and pass through the 
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Figure 7 Impact of the channel length scaling on the transfer characteristic for V GS = 0.5 V. 
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channel towards the drain, and then at that point when the 
charge goes to zero, the electrons are infused into the 
space charge region where they are swept by the E-field to 
the drain contact. Compared to the original length L, the 
change in channel length AL is small, then the drain 
current will be regular for Vd S > V DS (sat). The region of 
the Ir>-Vvs characteristic is referred to as the saturation 
region. When V G s is changed, the /d-Vds curve will also be 
changed. It was found that if V GS increases, the initial slope 
of IrrVns will be raised. We can also infer from Equation 14 
that the value of Vd S (sat) is a function of V G s- A family of 
curves is created for this n-channel enhancement-mode 
TGN SB FET, as shown in Figure 6. 

Also, it can be seen that by increasing V GS , the saturation 
current increases, showing the fact that a larger voltage 
drop occurs between the gate and the source contact. Also, 
there is a bigger energy value for carrier injection from the 
source contact channel [20]. The impact of power supply 
up-scaling decreases the SB length at the drain side, allow- 
ing it to be more transparent and resulting in more turn- 



on current to flow. Therefore, an acceptable performance 
comparable to the conventional behavior of a Schottky 
transistor is obtained. The scaling of the channel length 
improves gate electrostatic control, creating larger trans- 
conductance and smaller subthreshold swings. The effect 
of the channel length scaling on the I-V characteristic of 
TGN SB FET is investigated in Figure 7. It shows a similar 
trend when the gate-source voltage is changed. It can be 
seen that the drain current rises substantially as the length 
of the channel is increased from 5 to 50 nm. 

To get a greater insight into the effect of increasing 
channel length on the increment of the drain current, 
two important factors, which are the transparency of SB 
and the extension of the energy window for carrier con- 
centration, play a significant role [49,50]. For the first 
parameter, as the SB height and tunneling current are 
affected significantly by the charges close to the source 
of SB FET, the channel length effect on the drain current 
through the SB contact is taken into account in our pro- 
posed model. Moreover, when the center of the channel 





v Ds<voit> v Ds( voit) 



Figure 8 Comparison between proposed model and typical l-V characteristics of other types of transistors, (a) MOSFET with Si0 2 gate 
insulator [50] (V GS = 0.5V), (b) TGN MOSFET with an ionic liquid gate, C ins » C q [49] (V GS = 0.5 V), (c) TGN MOSFET with a 3-nm Zr0 2 wrap 
around gate, C ins ~ C q [49] (V GS = 0.37 V), (d) TGN MOSFET with a 3-nm Zr0 2 wrap around gate, C ins ~ C q [49] (V GS = 0.38 V). 
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Figure 9 Subthreshold regime of TGN SB FET at different values of V DS (V) for L = 25 nm. 
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of the SB FET is unoccupied with the charge impurities, 
the drain-source current increases because of the fact 
that free electrons are not affected by positive charges 
[49]. The effect of the latter parameter appears at the be- 
ginning of the channel where the barrier potential 
decreases as a result of low charge density near the 
source. This phenomenon leads to widening the energy 
window and ease of electron flow from the source to the 
channel [50]. Furthermore, due to the long mean free 
path of GNR [52-55], the scattering effect is not domin- 
ant; therefore, increasing the channel length will result 
in a larger drain current. 



For a channel length of 5 nm, direct tunneling from 
the source to drain results in a larger leakage current, 
and the gate voltage may rarely adjust the current. The 
transistor is too permeable to have a considerable dis- 
parity among on-off states. For a channel length of 10 
nm, the drain current has improved to about 1.3 mA. 
The rise in the drain current is found to be more signifi- 
cant for channel lengths higher than 20 nm. That is, by 
increasing the channel length, there is a dramatic rise in 
the initial slope of 7 D versus Vd S . Also, based on the sub- 
threshold slope model and the following simulated 
results, a faster device with opposite subthreshold slope 




0.5 0.6 
V GS (Voltage) 

Figure 10 / D (uA)-y GS (V) characteristic of TGN SB FET at different values of V DS . 
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or high on/off current ratio is expected. In other words, 
it can be concluded that there is a fast transient between 
on-off states. Increasing the channel length to 50 nm 
resulted in the drain current to increase by about 6.6 
mA. The operation of the state-of-the-art short-channel 
TGN SB FET is found to be near the ballistic limit. In- 
creasing further the channel length hardly changes nei- 
ther the on-current or off-current nor the on/off current 
ratio. However, for a conventional metal-oxide-semicon- 
ductor field-effect transistor (MOSFET), raising the channel 
length may result in the channel resistance to proportion- 
ally increase. Therefore, in this case, down-scaling the 
channel length will result in significant loss of the on/off 
current ratio as compared to the SG device. 

Figure 8 shows a comparative study of the presented 
model and the typical I-V characteristics of other types 
of transistors [49,50]. As depicted in Figure 8, the pro- 
posed model has a larger drain current than those tran- 
sistors for some value of the drain-source voltages. The 
resultant characteristics of the presented model shown 
in Figure 8 are in close agreement with published results 
[49,50]. In Figure 8, DG geometry is assumed for the 
simulations instead of the SG geometry type. 

In order to have a deep quantitative understanding of 
experiments involving GNR FETs, the proposed model is 
intended to aid in the design of such devices. The Si0 2 
gate insulator is 1.5 nm thick with a relative dielectric 
constant K = 3.9 [50] (Figure 8a). Furthermore, the gate- 
to-channel capacitance C g is a serial arrangement of in- 
sulator capacitance C ins and quantum capacitance C q 
(equivalent to the semiconductor capacitance in conven- 
tional MOSFETs). Figure 8b shows a comparative study 
of the presented model and the typical I-V characteristic 
of a TGN MOSFET with an ionic liquid gate. The avail- 
ability of the ionic liquid gating [49] that can be modeled as 
a wrap-around gate of a corresponding oxide thickness of 1 
nm and a dielectric constant e r = 80 results in C ins » C q , 
and MOSFETs function close to the quantum capacitance 
limit, i.e., C g ~ C q [49]. As depicted in Figure 8c,d, the com- 
parison study of the proposed model with a TGN MOSFET 
with a 3-nm Zr0 2 wrap-around gate for two different 
values of V GS is notable. A 3-nm Zr0 2 (e r = 25) wrap- 
around gate has C ins comparable to C q for solid-state 
high-A: gating, and this is an intermediate regime among 
the MOSFET limit and C q limit. 

Recently, a performance comparison between the 
GNR SB FETs and the MOSFET-like-doped source- 
drain contacts has been carried out using self-consistent 
atomistic simulations [20,21,48-50,56,57]. The MOSFET 
demonstrates improved performance in terms of bigger 
on-current, larger on/off current ratio, larger cutoff fre- 
quency, smaller intrinsic delay, and better saturation be- 
havior [21,50]. Disorders such as edge roughness, lattice 
vacancies, and ionized impurities have an important 



effect on device performance and unpredictability. This 
is because the sensitivity to channel atomistic structure 
and electrostatic environment is strong [50]. However, 
the intrinsic switching speed of the GNR SB FET is sev- 
eral times faster than that of the Si MOSFETs. This 
could lead to promising high-speed electronics applica- 
tions, where the large leakage of the GNR SB FET is of 
fewer concerns [20]. An efficient functionality of the 
transistor with a doped nanoribbon has been noticed in 
terms of on/off current ratio, intrinsic switching delay, 
and intrinsic cutoff frequency [48]. 

Based on the presented model, comparable with the 
other experimental and analytical models, the on-state 
current of the MOSFET-like GNR FET is 1 order of 
magnitude higher than that of the TGN SB FET. This is 
because the gate voltage ahead of the source-channel flat 
band condition modulates both the thermal and tunnel 
components in the on-state of MOSFET-like GNR FET, 
while it modulates the tunnel barrier only of the metal 
Schottky-contact TGN FET that limits the on-state 
current. Furthermore, TGN SB FET device performance 
can be affected by interlayer coupling, which can be 
decreased by raising the interlayer distance or mismatch- 
ing the A-B stacking of the graphene layers. 

It is also noteworthy that MOSFETs operate in the re- 
gion of subthreshold (weak inversion) as the magnitude 
of Vqs is smaller than that of the threshold voltage. In 
the weak inversion mode, the subthreshold leakage 
current is principally as a result of carriers' diffusion 
[58,59]. The off-state current of the transistor (/off) is 
the drain current when Vgs = 0. The off-state current is 
affected by some parameters such as channel length, 
channel width, depletion width of the channel, gate 
oxide thickness, threshold voltage, channel-source dop- 
ing profiles, drain-source junction depths, supply volt- 
age, and junction temperature [59]. 

Short-channel effects are defined as the results of scal- 
ing the channel length on the subthreshold leakage 
current and threshold voltage. The threshold voltage is 
decreased by reducing the channel length and drain- 
source voltage [58-61]. In the subthreshold region, the 
gate voltage is approximately linear [58,59]. It has been 
studied that the decrease of channel length and drain- 
source voltage results in shifting the characteristics to 
the left, and it is obvious that as the channel length gets 
less than 10 nm, the subthreshold current increases dra- 
matically [62]. Based on the International Technology 

Table 1 Subthreshold slope of TGN SB FET at different 



values of V DS 



V DS (mV) 


1 1.1 


1.2 1.3 


1.4 


1.5 


Subthreshold 


59.5238 54.1419 


49.6032 45.8085 


42.5134 


39.2542 


slope 










(mV/decade) 
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Roadmap for Semiconductors (ITRS) near-term guide- 
line for low-standby-power technology, the value of the 
threshold voltage is close to 0.3 V [59]. Figure 9 illus- 
trates the subthreshold regime of TGN SB FET at differ- 
ent values of drain-source voltage. As shown in this 
figure, for lower values of drain-source voltage, the 
threshold voltage is decreased and meets the guidelines 
of ITRS. 

The subthreshold slope, S (mV/decade), is evaluated 
by selecting two points in the subthreshold region of an 
^d-Vgs graph as the subthreshold leakage current is 
adjusted by a factor of 10. It has been noted that self- 
consistent electrostatics and the gate bias-dependent 
electronic structure have an essential role in determining 
the intrinsic limits of the subthreshold slope in a TGN 
SB FET, which stays well over the Boltzmann limit of 
the ideal value of 60 mV/decade or less than 85mV/dec- 
ade [58,63] .The subthreshold slope, as one of the key 
issues of deep-submicrometer devices, is defined as [59] 

y t - y off 

Subthreshold slope = , (15) 

log log /off 

where V t is the threshold voltage, V Q & is the off voltage 
of the device, /^ is the drain current at threshold, and 
/ 0 ff is the current at which the device is off. In other 
words, the subthreshold slope delineates the inverse 
slope of the log (/ D ) versus Vgs plotted graph as illu- 
strated in Figure 10. 

Average subthreshold swing is a fundamental param- 
eter that influences the performance of the device as a 
switch. According to Figure 10, the subthreshold slope 
for (/ = 100 nm) is obtained as shown in Table 1. 

Based on data from [64], for the effective channel 
lengths down to 100 nm, the calculated and simulated 
subthreshold slope values are near to the classical value 
of approximately 60 mV/decade. The subthreshold slope 
can be enhanced by decreasing the value of the buried 
oxide capacitance C BOX or by increasing the value of the 
gate oxide capacitance C GOX [64]. Based on the simu- 
lated results, it can be concluded that when the channel 
material is replaced by TGN, the subthreshold swing 
improves further. 

The comparison study between the presented model 
with data from [62,64] showed that due to the quantum 
confinement effect [39,43], the value of the subthreshold 
slope in the case of TGN SB FET is less than those of 
DG metal oxide semiconductor and vertical silicon-on- 
nothing FETs [62,64] for some values of drain-source 
voltage. A nanoelectronic device characterized by a steep 
subthreshold slope displays a faster transient between 
on-off states. A small value of S denotes a small change 
in the input bias which can modulate the output current 
and thus leads to less power consumption. In other 



words, a transistor can be used as a high-speed switch 
when the value of S is small. As a result, the proposed 
model can be applied as a useful tool to optimize the 
TGN SB FET-based device performance. It showed that 
the shortening of the top gate may lead to a considerable 
modification of the TGN SB FET current-voltage prop- 
erties. In fact, it also paves a path for future design of 
the TGN SB devices. 

Conclusions 

TGN with different stacking arrangements is used as metal 
and semiconductor contacts in a Schottky transistor junc- 
tion. The ABA-stacked TGN in the presence of an external 
electric field is also considered. Based on this configur- 
ation, an analytical model of junction current-voltage 
characteristic of TGN SB FET is presented. The depend- 
ence of the drain current versus the drain-source voltage 
of TGN SB FET as well as the back-gate and top-gate vol- 
tages for different values of gate-source voltage and geo- 
metric parameters such as channel length are calculated. 
In particular, we conclude that by increasing the applied 
voltage and also channel length, the drain current in- 
creases, which showed better performance in comparison 
with the typical behavior of other kinds of transistors. Fi- 
nally, a comparative study of the presented model with 
MOSFET with a Si0 2 gate insulator, a TGN MOSFET with 
an ionic liquid gate, and a TGN MOSFET with a Zr0 2 
wrap-around gate was presented. The proposed model is 
also characterized by a steep subthreshold slope, which 
clearly gives an illustration of the fact that the TGN SB 
FET shows a better performance in terms of transient be- 
tween off-on states. The obtained results showed that due 
to the superior electrical properties of TGN such as high 
mobility, quantum transport, ID behaviors, and easy fabri- 
cation, the suggested model can give better performance as 
a high-speed switch with a low value of subthreshold slope. 
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